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ABSTRACT 

We report the discovery of the Ukely IR counterpart to the Anomalous X-ray Pulsar (AXP) IRXS 
J 170849-4009 10, based on the combination of Chandra HRC-I X-ray position, and deep optical/IR observations 
carried out from ESO telescopes and the Canada France Hawaii Telescope (CFHT) during 1999-2002. Within the 
narrow uncertainty region we found two relatively faint (A''=20.0 and K'=n .SJ) IR objects. Based on their color 
and position in the J-K' versus J-H diagram only the brighter object is consistent with the known IR properties of 
the counterparts to other AXPs. No variability was detected for this source, similarly to what is observed in the 
case of 4U0142+614. Like in other AXPs, we found that the IR flux of IRXS J 170849-4009 10 is higher then 
expected for a simple blackbody component extrapolated from the X-ray data. If confirmed, this object would be 
the fourth IR counterpart to a source of the AXP class, and would make the IR excess a likely new characteristic 
of AXPs. 

Subject headings: stars: neutron — stars: pulsars: general — pulsar: individual: — IRXS J 170849-4009 10 — 
infrared: stars — X-rays: stars 



1. INTRODUCTION 

AXPs are thought to be solitary magnetic rotating neutron 
stars (NSs) either with a standard magnetised field or larger 
than 10'"* Gauss (magnetars), although the binary system sce- 
nario with a very low mass companion is not completely ruled 
out by current observational data, (for a review see Israel et 
al. 2002a; Mereghetti et al. 2002 and references therein). Dif- 
ferent production mechanisms for the observed X-ray emission 
have been proposed, involving either accretion or dissipation of 
magnetic energy. The magnetar model, originally proposed by 
Duncan & Thompson (1992) to explain Soft Gamma-ray Re- 
peaters (SGRs), appears to be successful at interpreting most 
of the properties of AXPs. In fact, AXPs have been linked 
to SGRs (thought to be magnetars, neutron stars powered by 
their strong magnetic fields) because of similar timing prop- 
erties, namely periods and positive period derivatives (Duncan 
& Thompson 1992; Thompson & Duncan 1993; 1996). What 
differentiates the two classes of objects is unclear The simi- 
larity in the spin parameters would not be sufficient by itself 
to differentiate AXPs and SGRs from other groups of pulsars 
with very different emission properties. Conversely, a very high 
magnetic field strength (if at all) cannot be the sole factor gov- 



erning whether or not a neutron star is a magnetar, a radio pulsar 
or in a binary system (Camilo et al. 2000). 

The possible connection of AXPs with SGRs has gained 
more credibility with the recent detection of SGR-like X-ray 
bursts from two AXPs (Kaspi & Gavriil 2002, Gavriil et al. 
2002) which also showed IR variabiHty. For IE 2259+586, IR 
variability of the likely counterpart has been detected few days 
after a strong bursting activity seen with the 7?XTE (Kaspi et 
al. 2002, Israel et al. 2003). Similarly, the variability of the 
IR counterpart to IE 1048.1-5937 is thought to be related to 
X-ray variability (Israel et al. 2002b), though no simultaneous 
X-ray/IR observations were available. The magnetar scenario 
does not make any prediction (yet) on the optical/IR emission 
seen in AXPs. If AXPs are powered by accretion via a disc, in 
analogy with the X-ray bursts observed in Low Mass X-ray Bi- 
naries, it is expected that the effects of the X-ray bursts activity 
in AXPs (and SGRs) propagate toward longer wavelengths e.g., 
via reprocessing in the disc (see Middleditch & Nelson 1976; 
Lawrence et al. 1983). Evidence for flattening (or excess) of 
the flux distribution in the IR band (with respect to a simple 
blackbody component forced to be consistent both with the X- 
ray and optical data) has been reported in three AXPs, namely 



' The results reported in this Letter are partially based on observations carried out at ESO, La Silla, Chile (63.H-0294, 66.D-0440, 67.D-0116 and 68.D-0350), and 
the Canada-France-Hawaii Telescope operated by the National Research Council of Canada, the Centre National de la Recherche Scientifique de France and the 
University of Hawaii (02bf21). 
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Fig. 1 .— Near-IR a:' band image of the region including tlie position of IRXS J17{)849-4()()9I(), taken from the 3.6-m CFHT. The ROSAT and Chandra uncer- 
tainty circles (9" and 0.8", respectively) are also shown. The inset shows the close-up of the Chandra position with the two counterpart candidates marked (letter 
A and B). Coordinates are RA (hms) and Dec. "; equinox J2000). 



1E2259+586, IE 1048.1-5937 and 4U0142+614 (Hulleman et 
al. 2001; Wang & Chakrabarty 2002; Israel et al. 2003). The 
spectral flattening and variability of the IR emission of AXPs 
are a potentially important diagnostics for the study of these 
enigmatic objects and their possible connections with other 
classes of pulsars. 

IRXS J170849-400910 was discovered by the ROSAT mis- 
sion during the All Sky Survey program (Voges et al. 1996), 
but only a few years later ~11 s pulsations were discovered 
by ASCA (Sugizaki et al. 1997). Both the period derivative 
(P~2xl0~^' ss~') and the optical upper limits inferred for its 
counterpart were found to be consistent with those of AXPs 
(Israel et al. 1999). The source is radio quiet (no pulsations de- 
tected in the radio band; Israel et al. 2002a) and a stable rotator 
comparable to radio pulsars (Gavriil & Kaspi 2002). It is also 
the only AXP for which X-ray spectral variations as a function 
of the pulse phase have been detected associated to variations 
of the pulse profile (similar to those observed in "standard" 
accretion-powered X-ray binary systems; Israel et al. 2003). 
Finally, IRXS J 170849-4009 10 represents the only AXP for 
which a spectral signature, thought to be due to resonant cy- 
clotron absorption, has been detected confirming the presence 
of a relatively high magnetic field (Rea et al. 2003). 

In this paper, we present the results obtained from Chandra, 
ESO and CFHT observations of IRXS J 170849-4009 10 in the 
X-ray, optical and IR bands. We identified the Ukely IR coun- 
terpart to this AXP based on positional coincidence and unusual 
IR colors. The overall (from X-ray to IR) energy distribution 
is also presented and discussed in the light of the IR excess de- 
tected in IRXS J170849-400910 and other AXPs. 

2. Chandra OBSERVATIONS 

The field of IRXS J 170849-4009 10 was observed during 
Chandra Cycle 2 with the High Resolution Imager (HRC-I; 



Zombeck et al. 1995) on 2001 September 23 for an effective 
exposure time of 9870 s. Data were reduced with CIAO version 
2.2 and analysed with standard software packages for X-ray 
data (CIAO, XIMAGE, XRONOS, etc.). The observation was 
carried out with a nominal aspect solution and the latest caU- 
bration files were used. Only one source was detected in the 
HRC-I (see Israel et al. 2002b for details on the detection al- 
gorithms). The source has the following coordinates: R.A. = 
17''08'"46?87, DEC. = ^0°08'52."44 (equinox J2000), with an 
uncertainty circle radius of 0."7 (90% confidence level; consis- 
tent with ROSAT HRI positions in Israel et al. 1999). Photon 
arrival times were extracted from a circular region with a ra- 
dius of l."5, including more than 90% of the source photons, 
and corrected to the barycenter of the solar system. A coher- 
ent pulsation at a period of 1 1 .00 1 1 ±0.0005 s (90% confidence 
level) was detected confirming that the source is indeed IRXS 
J170849-400910. 

We also carried out a spatial profile analysis for the (back- 
ground subtracted) X-ray emission from a region of 7" radius 
centered on the above reported source position. The spatial pro- 
file was found to be in good agreement with the expected Chan- 
dra Point Spread Function (PSF) for an on-axis source (see Is- 
rael et al. 2002b for details). Unfortunately, the high X-ray 
background level detected in the IRXS J 170849-4009 10 field 
prevented a more sensitive search for diffuse emission at larger 
radii. 

3. OPTICAL/IR OBSERVATIONS 

Most of the optical/IR data were obtained from the 3.6- 

m ESO telescope (ESO360; La Silla, Chile) equipped with 
ESO Faint Object Spectrograph and Camera 
(EFOSC2; a CCD with 2048x2048 pixels; 0."157 pixel scale 
and 5.2' X 5.2' field of view, FOV) in the optical band, and with 
the 3.5-m New Technology Telescope (NTT; La SiUa, Chile) 
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Fig. 2. — Color-color diagram obtained for the objects detected within a 
radius of 30" around the Chandra position of IRXS J 170849-4009 10. The 
proposed counterpart is marked (A). Star sequences are shown for compari- 
son, for different values of A\: 0, 2, 4, 6, and 11. 2 (thin solid Lines; the latter 
value obtained for the total Galactic absorption), and Av=7.8 (thick solid line; 
assuming A'h=1.4x 10^^ cm~^ as inferred from X-ray observations). The rect- 
angle which includes the position of object A represents the most probable 
region of the diagram where the IR counterpart to IRXS J170849^00910 
would lie, based on the known IR counterparts to AXPs. 

equipped with the Son OF Isaac (SOFI; Hawaii HgCdTe 
1024x1024 array; 0."292 pixel scale and 4.9'x4.9' FOV) in 
the near-IR band. Additional IR data were obtained at the 3.6- 
m CFHT (Mauna Kea, Hawaii) equipped with the Adaptive 
Optics Bonnette (AOB; Hawaii HgCdTe 1024x1024 ar- 
ray; 0."035 pixel scale and 36"x36" FOV). 

Optical Johnson R filter deep images were obtained with 
EF0SC2 on 1999 September 15-16 with effective exposure 
times of 1500s (seeing of 0"8). Standard reduction packages 
were used in the analysis of the optical data (DAOPHOT II; 
Stetson 1987) in order to obtain the photometry of each stel- 
lar object in the images. A limiting magnitude (3cr confidence 
level) of /?~26.5 was reached. 

Images in the / and H bands were initially acquired with 
SOFI on 2001 May 26 with 1920s {J) and 2400s (H) of ef- 
fective exposure time (seeing of 0"4-0"6). Observations yield 
a limiting magnitude of 22.7 (J) and 22.6 {H; 3(t confidence 
level). SOFI Ks images were obtained on 2002 February 19 
(seeing 0"8). The exposure time was of 3000 s, and a limiting 
magnitude of Ks^2\A was reached. In all these cases, single 
5 s-long exposure images were taken for each filter with offsets 
of 40" in order to sample and subtract the variable IR back- 
ground. Finally, H and K' images were obtained with AOB on 
2002 August 17 with effective exposure time of 2700 s (for each 
filter; seeing of 0."4-0."5). Thanks to the adaptive optics we ob- 
tained a source PSF of ~0."12 yielding a limiting magnitude of 
23.1 {H) and 21.8 {K'). During the latter run single 60 s-long 
exposure images were taken with random offsets in the 4"- 10" 
range. Standard IR software packages were used for sky frame 
subtraction and image coaddition (Eclipse and IRDR; Dev- 
illard 1997 and Sabbey et al. 2001). 

To register the Chandra coordinates of IRXS 
J 170849-4009 10 on our optical/IR images, we computed the 
image astrometry using, as a reference, the positions of stars 
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Fig. 3.— Broad band energy spectrum of IRXS J 170849-4009 10. X-ray 
raw data are taken from the medium and low energy instruments on board 
BeppoSAX (filled squares) while the solid upper curves are the unabsorbed 
fluxes for the black body (BB), the power law (PL), and the sum of the two 
components. On the lower left corner of the plot are the optical/IR fluxes: plus 
signs are absorbed values, while filled triangles and circles are unabsorbed 
ones for Av=7.8 and 11.2, respectively. The thick dotted, dash-dotted and 
solid lines are the fossil disk models discussed in the text. 

(about 20 for ESO instruments) selected from the GSC2.2 cat- 
alog which has an intrinsic absolute accuracy of about Of'2- 
0"4 (depending on magnitude and sky position of stars). Af- 
ter taking into account the uncertainties of the source X-ray 
coordinates (0"7), the rms error of our astrometry (0."12) and 
the propagation of the intrinsic absolute uncertainties on the 
GSC2.2 coordinates (we assumed a value of 0"3), we estimated 
the final accuracy to be attached to the IRXS J 170849-4009 10 
position of about 0."8. Figure 1 shows a region of 28"x20" 
around the IRXS J170849-400910 position in the K' band 
CFHT image (90% confidence level ROSAT and Chandra un- 
certainty circles superimposed). In the inset we show the close- 
up of the region around the Chandra position with the counter- 
part candidates marked. 

The 3.6-m ESO optical R band data show that no ob- 
ject is consistent with the Chandra position of IRXS 
J 170849-4009 10 down to a hmiting magnitude of about 26.5. 
On the other hand, two relatively faint objects are detected in 
the IR images carried out at the CFHT (marked with A and B 
in Figure 1). The latter two sources have the following IR mag- 
nitudes as inferred from CFHT (NTT) data: /:'=17.53±0.02, 
i/=18.85±0.05 {Ks=n3±Q.\, //=18.6±0.1 and 7=20.9±0.1) 
for object A, /s:'=20.00±0.08 and //=20.43±0.07 for B. Unfor- 
tunately, no J measurement is available for object B due to its 
faintness and likely confusion with object A in SOFT images. 

Figure 2 shows the color-color diagram of the objects in the 
field with star A marked. Note that the hydrogen column den- 
sity to IRXS J170849-400910, inferred from the X-ray obser- 
vations is A?H=l-40(6)x 10^^ cmr^ (Israel et al. 2001, Rea et al. 
2003), converts into a visual extinction of Ay ~ 7 .8 for a typical 
dust-to-gas ratio (Predehl & Schmitt 1995, see thick soUd line 
in Figure 2). In the J-K' versus J-H plane the reddening vec- 
tor Ues roughly parallel to the locus of points belonging to star 
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theoretical isochrones at solar chemical composition and stellar 
ages ranging from lOMyr to 8 Gyr (Bono at al. 1997, 2000; a 
value of 3.1 has been assumed for the parameter R=Av/E{B-Vy, 
Cardelli et al. 1989, Fitzpatrick 1999). We note that the syn- 
thetic sequences account for nearly all the objects but few red 
stars (J-K' colors larger than 3) that are uniformly distributed 
over the FOV. This is likely due to the parameter R, the value 
of which is Ukely larger than that assumed as the "nominal" 
one (in particular when looking at several spiral arms through 
the Galactic plane, as in the case of IRXS J 170849-4009 10). 
Finally, we used the information on the colors of the known 
IR counterparts to AXPs iH-K'=l.l and 1.4 for 4U0142+614 
and IE 1048.1-5937, respectively; Israel et al. 2003; Wang & 
Chakrabarty 2002), to identify the most probable region of the 
diagram where the IR counterpart to IRXS J 170849-4009 10 
would lie, if similar colors apply and assuming Ay >6.5. Object 
A is clearly included in the above region and it is also the only 
object showing unusual colors with respect to both the theo- 
retical star sequences and all other detected objects. All the 
above results make object A a robust candidate, and in the fol- 
lowing we will refer to it as the likely counterpart to IRXS 
J170849-400910. Based again on the color, we regard candi- 
date B as a far less likely candidate {H-K'=OA). 

4. DISCUSSION 

The proposed IR counterpart to IRXS J170849-400910is by 

far the brightest one among those of AXPs: K'=21.7, K'=20.0, 
/i:.9=19.4 for 1E2259+586,4U0142+614, and IE 1048.1-5937, 
respectively (Hulleman et al. 2001 and 2002, Israel et al. 
2003, Wang & Chakrabarty 2002). We note that the rela- 
tively large X-ray-to-IR unabsorbed flux ratio of about 500 
is lower than that obtained for other AXPs. However, IRXS 
J170849-400910 is also the most luminous AXP in the 0.5- 
lOkeV band with an unabsorbed luminosity of > 4x10^^ 
ergs"^ (assuming a distance >5kpc). In order to characterise 
the broadband energy distribution of the source, we plotted the 
IR through X-ray data of IRXS J 170849-4009 10 in Figure 3. 
X-ray fluxes have been inferred by using the phase-averaged 
spectral parameters obtained with BeppoSAX data (Israel et al. 
2001, Rea et al. 2003). Two different values of the extinction 
have been used for model fitting purposes: Av=7.8 and 11.2 
(triangles and circles, respectively), corresponding to the Nh 
inferred from the X-ray spectra and the Galactic absorption in 
the direction of the source, respectively. 

The blackbody component detected in the X-rays (kT ^0.45 
keV; Israel et al. 2001, Rea et al. 2003) cannot account for 
the (relatively high) IR flux. Similarly, the power-law compo- 
nent, if extrapolated to IR wavelengths, would imply a much 
higher IR flux level which is simply not observed; a cut-off 
must therefore occur somewhere in the UV/optical bands. Re- 
gardless of the exact position of the cut-off and its origin, we 
note fliat the X-ray-to-IR emission of IRXS J 170849-4009 10 
cannot be fit with a simple spectral component (similar to the 
case of 4U 0142-1-614 and IE 2259-1-586; Hulleman et al. 2000, 
2001). 

The magnetar model does not account for the observed IR 
emission or IR variabiUty of AXPs and no predictions can be 
verified. Therefore, we discuss the above observational findings 
in the context of models based on fossil disks (Chatterjee et al. 
2000; Alpar 2001) since they make clear predictions for the IR 
emission. However, up to now the latter models are unable to 
account for the bursts seen in AXPs and SGRs.. In the case of 



a fossil disk the IR emission arises mainly from two compo- 
nents: viscous dissipation and reprocessing of X-ray flux im- 
pinging on the disk from the pulsar (Perna et al. 2000, Perna & 
Hernquist 2000, Alpar 2001). Specifically, at IR wavelengths 
and based on the X-ray luminosity of IRXS J 170849-4009 10, 
the dominant emission component is expected to be X-ray irra- 
diation. The models we considered have been computed for a 
distance of 5 kpc and the observed X-ray luminosity as an input 
parameter, but since the disk luminosity scales almost linearly 
with the X-ray luminosity, the result is nearly independent of 
the distance. As far as the X-ray spectrum, in the accretion 
model (with magnetic fields of ~10'^ Gauss) one would expect 
the emission to be roughly a blackbody produced in a region 
(polar cap) much smaller than the area of the star Current X- 
ray data are consistent with this, but not good enough to allow 
a firm discrimination with respect to the X-ray spectrum of a 
magnetar, in which case one would expect the area of the emit- 
ting region to be consistent with the whole surface of the star 
(Perna et al. 2001). 

In Figure 3 we show disk models obtained with the inner and 
outer radius of the disk, r,„ and rout, and the inclination left 
free to vary (however, note that the inclination only scales the 
fluxes, while r,„ and ro„, modify the shape of the model and the 
flux peak position; see Perna et al. 2000 for details). The best fit 
was obtained for a smaU truncated disk, with r,„ and rgut equal 
to 10" and 5xlO" cm, respectively (dotted Une in Figure 3). 
However, this solution would not be able to explain the X-ray 
emission from AXPs, being the value of r,„ larger than the coro- 
tation radius (~8 x 10^ cm), the characteristic value of r,„ in or- 
der to have accretion on the neutron star surface, and at the 
same time, spin-down. Therefore, we fixed r,„ to the corotation 
radius for two different values of extinction. The results of the 
fit are shown in Figure 3 (dashed-dotted and solid lines, respec- 
tively). Values of ?'oi,f=10'^ and 8x10'^ cm were obtained for 
Ay=7.8 and 11. 2, respectively, (fallback compact disks; Menou 
et al. 2001). In both cases the / measurement is lower than the 
model expectations (2a). The above solution might be applied 
also to the case of an accretion disk formed by mass transfer 
from a (light) companion, although in this case an (unknown) 
amount of IR flux should be originated also from the star sur- 
face. We note that IR flattening has been recently discovered 
also for 4U 0142+614, and likely during an X-ray bursting ac- 
tivity phase of IE 1048.1-5937 (Wang & Chakrabarty 2002, 
Israel et al. 2003). The IR flattening or excess might represent 
a new important property shared by AXPs. 

The IR emission of AXPs in the magnetar scenario, remains 
to be addressed through detailed modehng. However, we note 
that the presence of a disk would not be necessarily in con- 
trast with the magnetar scenario. Three AXPs are associated 
with supernova remnants and likely embedded in dense re- 
gions. Therefore a "hybrid" scenario in which a magnetar ir- 
radiates a fossil disk (the matter of which is not accreting in 
this case) might not be unreasonable, though it should be con- 
firmed through additional future observations and its consis- 
tency checked through theoretical studies. 

In conclusion, even though we propose a possible solution 
for the IR excess observed in AXPs, we note that none of 
the proposed theoretical models (at least in their present form) 
seem to be able to easily account simultaneously for the IR, 
optical and X-ray emission of AXPs. 
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